Ground state and edge excitations of quantum Hall liquid at filling factor 2/3 
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We present a numerical study of fractional quantum Hall liquid at Landau level filling factor v — 
2/3 in a microscopic model including long-range Coulomb interaction and edge confining potential, 
based on the disc geometry. We find the ground state is accurately described by the particle-hole 
conjugate of a = 1/3 Laughlin state. We also find there are two counter-propagating edge modes, 
and the velocity of the forward-propagating mode is larger than the backward-propagating mode. 
The velocities have opposite responses to the change of the background confinement potential. On 
the other hand changing the two-body Coulomb potential has qualitatively the same effect on the 
velocities; for example we find increasing layer thickness (which softens of the Coulomb interaction) 
reduces both the forward mode and the backward mode velocities. 



PACS numbers: 

I. INTRODUCTION 

The fractional quantum Hall (FQH) effect is a remark- 
able phenomenon observed in two-dimensional electron 
gases (2DEGs) in a strong perpendicular magnetic field. 
FQH liquids are gapped and believed to possess topolog- 
ical order. ^ In particular, it has been established that the 

= 1/3 Laughlin state represents an Abelian topologi- 
cal phase. The excitations in such a phase can carry a 
fraction of electron charge and have fractional statistics 
which are in-between bosonic and fermionic statistics. 
Various experiments have reported observation of frac- 
tional charge^i^ii^. Recently, a series of experiments^'^ ob- 
served the so-called superperiods in the conductance os- 
cillations in an FQH quasiparticle interferometer, which 
have been interpreted as a reflection of fractional statis- 
ticsi^ii^ The bulk topological order is also reflected in the 
corresponding edge excitations, which are gapless. For 
a u = 1/S FQH liquid, with a sharp confining potential 
(no edge reconstruction) , there is only a single branch of 
bosonic excitations at the edge. The bosonic edge mode 
is chiral, i.e., propagating along the edge in one direc- 
tion (determined by the E x B drift) only, because the 
magnetic field breaks the time-reversal symmetry. The 
edge physics can be described by the chiral Luttinger liq- 
uid theory and has been verified by numerical tests in 
microscopic modelsj^i^^ 

In a hierarchical state, the FQH liquid supports mul- 
tiple branches of edge excitations. Depending on the 
bulk topological order, the edge modes may propagate 
in the same direction or in opposite directions. The sim- 
plest case with counter-propagating edge modes is the 
spin-polarized FQH liquid at filling fraction v = 2/3, 
which can be regarded as the particlc-hole conjugate of a 

= 1/3 Laughlin state or, equivalently, a hole Laughlin 
state embedded in a ;/ = 1 integer quantum Hall back- 



ground (see Fig. [T] for an illustration) Roughly 
speaking, the inner and outer edges are located at density 
changes of 2/3 1 and 1^0, respectively. In general, 
the two edge modes are coupled to each other and their 
properties may be dominated by disorder in the presence 
of random edge tunneling.— The edge physics of the 
I' = 2/3 state is intriguing since one of the edge modes 
propagates opposite to the classical skipping orbits dic- 
tated by the uniform magnetic field, leading to a neg- 
ative contribution to thermal Hall conductivity.^^ The 
counter-propagating edge modes have since been stud- 
jg(^ji8ji9j20j2^ both theoretically and experimentally in 
recent years. 

Recently, a similar but more delicate situation arises at 
filling fraction v = 5/2, where the Moore- Read Pfaffian 
stateS^ and its particle-hole conjugated state, dubbed 
the anti-Pfaffian statC j^'^i^^ compete for the ground state. 
In the absence of Landau level mixing, impurity, or 
edge confinement, the particle-hole symmetry is unbro- 
ken. In this case the two states, in the bulk, are ex- 
pected to be degenerate in the thermodynamic limit. 
But these two states have very different edge struc- 
tures: the Pfaffian state supports two co-propagating chi- 
ral edge modes (a charged bosonic mode and a neutral 
fermionic mode), while the anti-Pfaffian state supports 
three counter-propagating charge and neutral modes i^i^i 
Their relation is somewhat similar to that between the 
I' = 1/3 and v — 2/?> edge states. We note very re- 
cent experiments2^'^'21 found indications of quasiparti- 
cle excitations with charge e/4 supported by both states 
and, interestingly, tunneling experiments^^ seem to favor 
the anti-Pfaffian state over the Pfaffian state. We also 
note that the particle-hole conjugates of the Read-Rezayi 
state^^ have been studied theoretically, with emphasis on 
properties of their edge excitations^^ 

Motivated by the recent work on the particle-hole con- 
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jugate of the Pfaffian state^^i^i and Read-Rezayi states?^ 
as well as by the experimental measurement of the I — V 
spectroscopy between individual edge channels,— we re- 
visit the polarized z/ = 2/3 FQH state with a detailed 
numerical study on the edge modes of the v — 2/3 
FQH droplet using a semi-realistic microscopic model. 
We find the ground states of the system for a wide pa- 
rameter range are well described by the composite of a 
h' — 1 IQH droplet and a z/ = 1/3 Laughlin hole droplet. 
The number of electrons and holes in the two droplets 
vary as the strength of confining potential varies, under 
the constraint that the total number of electrons does 
not change. Two counter-propagating edge modes are 
clearly visible in our results. Quantitatively, we find 
the forward-propagating outer edge mode (arising from 
the IQH edge) has a larger velocity than that of the 
backward-propagating inner edge (from the hole FQH 
edge). The structure of the excitation spectrum of the 
inner edge is identical to that of the Laughlin state at 
I' = 1/3 except for direction of propagation. Increas- 
ing the edge confining potential increases the outer edge 
mode velocity and reduces the inner edge mode veloc- 
ity. We also carry out a particle-hole transformation of 
the electronic Hamiltonian with hard-core interaction to 
generate the Hamiltonian that makes the hole Laughlin 
state and hole edge states as its exact zero-energy ground 
states. Using a mixed Hamiltonian which contains both 
the electron Coulomb Hamiltonian and the the conju- 
gate Hamiltonian of the two-body hard-core interaction, 
the bulk excitation energies can be raised to allow for 
a clearer separation between bulk and edge excitations. 
We find our results are robust in the presence of the elec- 
tronic layer thickness, whose main effect is softening the 
Coulomb interaction and reducing the velocities of both 
edge modes. 

The rest of the paper is organized as follows. In Scc.HIl 
we describe the microscopic model used in this work. We 
discuss the nature of the ground states in Sec. IIIII We 
present the overlap study of the ground states with vari- 
ational wave functions in Sec. IIVI A detailed analysis 
of the edge states follows in Sec. |Vl One can find ex- 
plicit construction of the Hamiltonian for the variational 
wave functions of the ground states and edge states in 
Sec. I VII We consider the effect of electron layer thick- 
ness in Sec. I VIII and quasihole excitations in Sec. IVIIII 
Finally, some concluding remarks are offered in Sec. IIXI 



II. THE MODEL 

In recent years, we have developed a semi-realistic mi- 
croscopic model for FQH liquids and have studied edge 
excitations and instabilities, quasiholc/quasiparticle ex- 
citations, and edge tunneling in Laughlin and Moore- 
Read Pfaffian phases .iSiSi^iSS^^ rpj^^ advantage of 

the model is that, depending on the parameters, the 
Laughlin phase, the Moore-Read phase, as well as edge 
reconstructed states and quasihole/quasiparticle states, 



emerge naturally as the global ground state of the mi- 
croscopic Hamiltonian without any explicit assumptions, 
e.g., on the value of the ground state angular momentum. 
This way, we can study the stability of phases and their 
competition. Another advantage of the model is that we 
can analyze the edge excitations of the semi-realistic sys- 
tem and identify them in a one-to-one correspondence 
with edge excitations of the corresponding edge theory 
(or conformal field theory). In addition to confirming 
the bulk topological order, we can use the microscopic 
calculation to extract energetic quantities, such as edge 
velocities, which are crucial for quantitative comparisons 
with experiments; for example in a recent study we used 
the edge mode velocities extracted from our numerical 
study to estimate the quasiparticle dephasing length at 
finite temperatures a,t v = 5/2.^^ In this paper, we apply 
our model and methods to the v = 2/3 FQH system. 

In this model, we consider a 2DEG confined to a plane 
with rotational symmetry. There is a neutralizing back- 
ground charge distributed uniformly on a parallel disk of 
radius i? at a distance d above the 2DEG (see Fig. 1 of 
Ref. [l3 for an illustration). The total charge of the disc is 
N^e, where is the number of electrons confined to the 
plane and the radius R or, equivalently, the density of the 
background charge is determined by the filling fraction 
v. We consider v = N^/N^ = 2Nel%/R'^, where Ib is the 
magnetic length and is the number of flux quanta 
enclosed in the disc. The distance d parameterizes the 
strength of the electron confining potential due to attrac- 
tion from background charge, which becomes weaker as 
d increases. We assume the electrons are spin-polarized, 
which is the case in strong magnetic fields. In the second 
quantization language, the Hamiltonian is written as: 



{mi} rn 



where the Coulomb matrix elements V{„.} are 



(1) 



Vlnu} " y ^ y ^ ?'20™i(ri)0*„2(^2)^^0„i3(ri)(/)„i4(f2), 

(2) 

and the background confining potential as a function 
of d is 



(3) 



Here e is the dielectric constant. We use the symmetric 
gauge A = ■^)] the single-particle wave function 

dm in the lowest Landau level is: 



0„,(z) = (27r2™m!)-i/2z™e-l^l'/'*. 



(4) 



Throughout the paper, we use the magnetic length Is as 
length unit and e^/eZs as energy unit. 

It is often convenient to cast the Coulomb matrix el- 
ements into a weighted sum of pseudopotentials intro- 
duced by Haldane.'^^ One of the advantages of express- 
ing two-body interactions in terms of pseudopotentials 
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FIG. 1: (color online) Schematic picture of electron density 
profile along radial direction of a = 2/3 FQH droplet. We 
assume there is a hole Laughlin state embedded in the v = 
1 electron background. Therefore, there are two interfaces: 
one is between the 1/3 hole Laughlin state and the integral 
filling background and the other is between the v = \ integer 
quantum Hall droplet and the vacuum. 

is that the Laughlin states become exact ground states 
for specific pseudopotential Hamiltonians; for example 
at filling fraction v — 1/3, with hard-core interaction be- 
tween electrons (in Haldane's pseudopotential language, 
= Si^m) and in the absence of confining potential, the 
Laughlin stated 

N ^ N 

*i/3(^i, ■■■ZN)^Y[{z^- z,f exp{-- ^ (5) 

is the exact ground state with zero energy, which exists in 
the subspace of total angular momentum M = 3Ne{Ne — 
l)/2, for Ne electrons in at least Norb = SiVg — 2 orbitals. 

III. GROUND STATE QUANTUM NUMBERS 

To begin our study on the v = 2/3 system, we ask 
to what extent we can conclude that the ground state of 
the semi-realistic model can be described by the particle- 
hole conjugate of the 1/3 Laughlin state on an IQH back- 
ground. The schematic profile of the electron density is 
shown in Fig.[l]in which we neglect the density oscillation 
of the Laughlin state for holes near its edge (see realistic 
curves in Fig. [4] for details). Suppose the system contains 
Ne electrons filling up to the Nj =th orbital (with single 
particle angular momentum m — 0,1, TV/— 1). Accord- 
ing to this picture, we have two droplets: Nj electrons 
fill the lowest Landau level (LLL) and form a lyj — 1 IQH 
state; in addition, Nh = {Nj — N^) holes form a, i^h — 1/3 
hole Laughlin state. The total angular momentum for 
such a state is 

^ Nj{Nj - 1) _ 3{Nj-N,){Nj^N,-l) 

2 2 ■ \ ) 

To reveal such a state in a microscopic calculation, one 
needs Norb > Nj orbitals. For example, with A^e = 20 
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FIG. 2: (color online) Phase diagram for systems with 20 
electrons at filling factor = 2/3, as a function of number of 
single electron orbitals Norb and background charge distance 
d. The ground state total angular momentum M changes 
as d increases. Angular momenta M — 270, 280, 288, and294 
correspond to 5-, 6-, 7-, and 8-hole Laughlin ground states 
respectively. However the M=294 state with Norb = 27 
is expected to be a stripe state which can be described by 
|*sp) = lOllllllOOOOOOllllllllllllU). 



electrons filling Nj — 26 orbitals, we have a 6-hole Laugh- 
lin droplet (filling the innermost 16 orbitals) and a vi = 1 
IQH droplet (filling all 26 orbitals). The total angular 
momentum of the state is M — 280. We note that for 
the specific system the average filling fraction is some- 
what different from i/ — 2/3, due to the higher density 
near the edge. 

If the hole-droplet picture is correct, one can possi- 
bly (but not necessarily for energetic reasons) find that 
the global ground state of the semi-realistic model has 
a total angular momentum M = 280, if 20 electrons are 
distributed in Norb > Nj = 26 orbitals. 

We study the global ground state of a system of 20 
electrons with various d and Norb (which serves as an 
additional hard edge confinement). We plot the results 
in Fig. O When the distance of the background charge 
d increases, the confining potential becomes weaker and 
the total angular momentum of the global ground state 
increases and goes through steps at M = 270, 280, 288, 
and 294. According to Eq. ([5]), these states correspond to 
Nh — 5, 6, 7, and 8 holes, respectively. Unlike the Laugh- 
lin or Moore- Read Pfaffian cases where M is uniquely de- 
termined by the number of electrons and a change of M 
is an indication of instabilityj^i^^ here we have a series 
of "good" ground states, since the number of electrons 
only fixes the difference between the IQHE droplet and 
the hole droplet. In other words, the two edges can move 
simultaneously to respond to the change of edge confin- 
ing potential. In this paper, we will focus on the 6-hole 
state, although the results should be general enough for 
other cases. 

The agreement between the actual series of ground 
state angular momenta and the prediction of Eq. ^ 
suggests that the picture of a Laughlin hole droplet em- 
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bedded in an IQH droplet is a good description of the 
ground states dX v — 2/3. In the foUowing sections we 
will present more evidence that this is indeed true for 
most cases. However, we would like to point out an ex- 
ception here. The quantum number M = 294 is con- 
sistent with an 8-hole state. But according to earlier 
discussion, we need Norb ^ 28 orbitals to accommodate 
this state. The reader may have already noticed that 
it also appears in the Norb — 27 case. This contra- 
diction suggests that there is another state (or phase) 
that is competing with the hole-droplet picture. In fact, 
it even suppresses the 7-hole state (M = 288). Analy- 
sis of the electron occupation number suggests that this 
state is likely a stripe state. We can use a fermionic 
occupation number configuration to represent the state 
as \^sp) = lOllllllOOOOOOllllllllllllll), where the 
string of Os and Is represents the occupation number of 
the corresponding single-particle orbital (from the left 
m = 0, 1, Norb — 1)- The stripe state is more compact 
than the 8-hole state and is expected to be energetically 
more favorable for the restricted case with Norb = 27 . 
To confirm this and to quantify the trend from Nort — 27 
to 30, we calculate the overlap between the ground state 
with M = 294 and l^'sp) and find that the ground state 
indeed has a large overlap (|(*sp|vE'm=294)P = 25.4%) 
with a stripe state when Norb — 27 and d — 1.21b- This 
overlap decreases as we increase the number of orbitals 
to 19.4% in 28 orbitals, 16.2% in 29 orbitals, and 8.2% in 
30 orbitals, suggesting the ground state for large enough 
Norb is not the stripe phase, but possibly the 8-hole state. 
Note that for well-defined overlap, we need to add the 
proper number of Os to the right of the occupation num- 
ber configuration |5'sp). 

IV. GROUND STATE WAVE FUNCTIONS AND 
OVERLAPS 

In this section we study the ground state wave func- 
tions at J/ = 2/3 and show that they arc consistent with 
the picture of a = 1/3 hole droplet on top of a i/g = 1 
electron droplet. For comparison between the two, we 
construct variational wave functions for u = 2/3 ground 
states by particle-hole conjugation of electron Laughlin 
states and calculate overlaps between them. To be spe- 
cific, we target the 6-hole ground state with M = 280. 
We first fill the lowest 26 single-particle LLL orbitals 
to obtain an IQH state, which can be represented by a 
string of single-particle occupation numbers |111---111). 
Then, we construct a 6-hole Laughlin wave function in 
the following manner. A 6-electron Laughlin state, which 
is partially occupying the lowest 16 orbitals, can be writ- 
ten as 

\Lte)^ n (^^-^.)'' (7) 

l<i<j<6 

where we have omitted a normalization constant and the 
Gaussian factor exp{— ^ • |zip/4}. We note that the 



Norb I 26 27 28 29 30 

Size of HS 1123 10867 54799 184717 473259 

\{^gs\Lle)fi\^J'^ 0.9401 0.7868 0.7012 0.6431 0.6011 

TABLE I: Overlaps between 20-electron ground states with 
M = 280 and d = 0.7Zb and the particle-hole conjugate of the 
6-electron Laughlin state obtained from the hard- core Hamil- 
tonian for several different numbers of orbitals Norb- In the 
case of d = 0.7, the ground state with Af = 280 is the global 
ground state for all cases in Fig. [21 The Hilbert subspace (HS) 
size for M = 280 increases rapidly when Norb increases from 
26 to 30 by about 450 times; however the overlap decreases 
slowly, indicating the robustness of the state. 

Laughlin state is a many-body wave function and can be 
written, symbolically in the second quantization form, as 

\L%) = E (8) 
Ki 

where 

|0)i6 = I000-- -000)16 (9) 

is the vacuum in 16 orbitals. Therefore, the 6-hole 
droplet embedded in 16 orbitals is 

l-^ie) = ^ a[i^T^Ci^Ci^Ci^Ci^Ci,^Ci^\ni ■■■ 111) 16, (10) 

which contains 10 electrons. After adding the additional 
10 filled orbitals, we have a many-body variational wave 
function for 20 electrons in 26 orbitals, which we denote 
as \LfQ)l° = iLfg) «) |llllllllll)io. For Norb > 26, we 
can add trailing Os for the empty orbitals at the edge 
accordingly, which we denote as l^ig)^ ^- We compare 
l^i6)wo^b to the ground state of our semi-realistic model 
with d = 0.71b, at which all ground states for Norb = 26- 
30 have M = 280, as shown in Fig. [H Table U shows the 
overlap between the global ground state \'^gs) in differ- 
ent numbers of orbitals and the 6-hole variational wave 
function [ifg)^ ^. As Norb varies from 26 to 30, the 
dimension of the Hilbert space increases by a factor of 
about 450, while the overlap still survives at about 60%. 
The decrease is largely due to the fact that the outer 
edge is no longer sharp as the angular momentum cut-off 
{Norb — 1) increases. 

In addition, the strength of the confining potential due 
to neutralizing background charge also affects, though in 
a minor way, the overlap between the v = 2/3 electron 
ground state and the variational wave function |-/jig)|J' ^. 
Fig. [3] shows that, for the rcginc in which the total angu- 
lar momentum of the global ground state is M — 280, the 
overlap for 20 electrons in 26 orbitals decreases from 0.94 
to 0.93 as d increases from 0.67 to 1.2. As the distance d 
between the 2DEG and the background charge increases, 
the confining potential the electrons experience becomes 
weaker. The electron wave function can expand, leading 
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FIG. 3: Overlap between the 20-electron ground state in 
26 orbitals and the particle-hole conjugate of the 6-electron 
Laughlin state as a function of d. The decrease of the over- 
lap indicates that the particle-hole conjugate Laughlin state 
favors smaller d, i.e., stronger confinement. 
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FIG. 4: (color online) Density profiles for the 6-electron 
Laughlin state, the 20-electron ground state in 30 orbitals 
with M = 280 and d = 0.71b, the sum of them, and the 
26-electron IQH state. The sum is almost the same as the 
density of IQH. 



to a smaller overlap. This is consistent with the Norb de- 
pendence, in the sense that the hole-droplet ground state 
favors stronger confinement. But since the range of d for 
a ground state with a certain number of holes is narrow, 
we can neglect the 1% change. 

An alternative way to compare states is to contrast the 
electron density profiles (though it is possible that two 
wave functions with identical density profiles can be or- 
thogonal to each other) . In Fig. [4l we plot the density 
profile of a 6-electron Laughlin state and that of the 6- 
hole ground state obtained by exact diagonalization of a 
system of 20 electrons in 30 orbitals at d = 0.71b- To 
compare, we plot the sum of the density profiles of the 
Laughlin state and the 6-hole state, together with the 
density profile of an IQH state with 26 electrons. It is 
clear that the density sum is almost the same as the den- 



sity of the IQH state, except for small oscillations which 
can be attributed to the long-range Coulomb interaction 
in the semi-realistic microscopic model. 

V. EDGE EXCITATIONS 

The analysis that can further confirm the picture that 
the V — 2/3 state is a Laughhn hole-droplet embed- 
ded in an IQH background is the study of edge exci- 
tations. In topological systems like FQH liquids, edge 
states have been demonstrated to be very effective and 
essential probes of the bulk topological order in both the- 
oretical calculations and experiments. As shown in FigdJ 
one expects two counter-propagating edge modes origi- 
nating from the two edges at 2/3 ^ 1 and 1 9 ,11,12,13 
In the same spirit as the analysis in Refs. [s^ and [13, we 
can label each low-energy edge excitation by two sets of 
occupation numbers {nL{lL)} and {nij(^ij)} for the inner 
and outer edge modes with angular momenta II, Ir and 
energies e^, e/j, respectively, nj^ili/) and nii{lfj) are non- 
negative integers. The angular momentum and energy of 
an edge excitation, measured relatively from those of the 
ground state, are 

AAf = - Y,^L{h)lL + Y.'^R^^R)^R^ (11) 
AE = ^nL(/L)eL(?i)+^nfl(/fl,)efl(/;i). (12) 

For the latter we assumed absence of interactions among 
the excitations. The negative sign in Eq. pT|) indicates 
the inner edge mode is propagating in the opposite di- 
rection to the outer edge mode. 

In Fig. O we plot the low-energy spectrum for 20 elec- 
trons in 28 orbitals sX d = 0.51b, whose global ground 
state has M = 280. The edge states are labeled by red 
bars based on the analysis we will discuss in the following 
paragraphs. Here we first point out that the inner edge 
excitations have negative AAf and are separated by an 
energy gap from other states (presumably bulk states) in 
each angular momentum subspace. The number of these 
inner edge states (including the ground state) are 1, 1, 
2, 3, and 5 for AM = 0-4, as predicted by the chiral 
boson edge theoryJ^ They have significant overlap with 
a 6-hole Laughlin droplet with corresponding edge exci- 
tations embedded in a 26-electron IQH background. On 
the other hand, the outer edge excitations {AM > 0) 
have higher excitation energies and are mixed with bulk 
states. In particular, for AM = 1, the edge state is the 
second lowest eigenstate in the M = 281 subspace. The 
state has a large overlap (62.7%) with the 6-hole Laugh- 
lin droplet embedded in the 26-electron IQH state with 
an edge excitation at AM — 1. Obviously, the outer edge 
mode has a larger velocity than the inner edge mode, con- 
sistent with the different charge density associated with 
the edge modes. 

To identify the inner edge excitations, we compare the 
edge spectrum to that of a 6-electron system at v = 1/3, 
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FIG. 5: (color online) Low-energy spectrum for 20 electrons in 
28 orbitals with d = 0.5. The inner and outer edge states are 
labeled by red bars and the simplest combination of the two 
is labeled by a black bar. Its energy (0.1049) is roughly the 
sum of the two lowest excitation energies for the two modes: 
0.02807(inner)+0.07687(outer). This state has a moderately 
large overlap(0.3959) with the particle-hole conjugate of a 6- 
electron Laughlin state with AAI = 1 edge excitation and 
embedded in the IQH edge state with AM = 1: iLf^iAM = 
1)---1101 >. 



with the neutralizing background charge at the same dis- 
tance from the 2DEG. In Fig. [HI we plot side-by-side 
the u = 1/3 and the i/ = 2/3 edge spectra, both with 
d — 0.51b- The edge states are labeled by red bars. 
From the comparison, one clearly sees that the u — 2/3 
state and the v = 1/3 state have similar edge excita- 
tions, but along opposite directions. The one-to-one cor- 
respondence of the edge excitations in the two cases can 
be established by studying overlaps of the correspond- 
ing pairs. Of course, the overlap that we really calcu- 
lated is the overlap between the eigenstate for — 2/3 
and the corresponding particle-hole conjugated state for 
1^ — 1/3 embedded in the 26-electron IQH background. 
To minimize the influence from the momentum cut-off, 
we choose 6 electrons in 22 orbitals with the same back- 
ground {d = 0.51b) for 1/3 filling. The results of the over- 
lap calculation are summarized in Table [ill The overlap 
becomes smaller as we go to higher energy, but remains 
above 40% up to |AAf| = 4. Overlaps between other 
pairs of states are significantly smaller. 

In order to identify edge excitations with positive AM, 
we need to consider the IQH edge excitations of the outer 
edge. A 26-electron IQH ground state in 28 orbitals 
can be represented by occupation numbers |11 • • • 11100), 
with 26 consecutive Is followed by two Os. So after 
adding a 6-hole Laughlin hole droplet, we denote the 
ground state as |Zfg • • • 11100) for convenience, although 
the many-body variational state cannot be written as 
a single occupation number string (i.e., a Slater deter- 
minant). In the same spirit, we can construct and de- 
note variational wave functions with IQH edge excita- 
tions as |Zfg • • • 11010) for AM = 1, \L% ■ ■ ■ 11001) and 



FIG. 6: (color online) Comparison of the low energy excitation 
spectrum of the u — 1/3, 6-electron Laughlin edge states (6 
electrons in 22 orbitals with Coulomb interaction, d = 0.51b) 
(a) and i/ = 2/3, 6-hole Laughlin edge states(20 electrons in 
28 orbitals, d = 0.51b) (b). The edge states are labeled by 
red bars. The overlaps between these edge states are shown 
in table nil 

jZfg • • • 10110) for AM = 2 in order to emphasize the ex- 
citation at the outer edge. We calculate the overlap be- 
tween the variational wave function |L^g • • • 11010) and 
the states in M = 281 subspace. The largest we find 
is the second state, |(*M=28i2|Z?6---11010)P = 0.6272, 
which we identify as the outer edge excitation. Simi- 
larly for AM = 2, we identify the fifth state, which has 
|(*A/=2825|i?6--10110)p = 0._1896, and the seventeenth 
state, which has |(^'a//=282i7 |ii6---11001) P = 0.1794, as 
edge excitations. We note the overlap already becomes 
small for AM = 2, indicating significant mixing between 
the edge states and bulk states. This is not surprising 
since, for the small system we consider, there is no gap 
protecting the edge states. 

Eqs. pT|) and p2)) suggest that there are also com- 
posite excitations that are combinations of these two 
counter-propagating edge modes. The simplest one is 
the combination of the edge states with AM — —1 and 
AM = 1, which resides in the M — 280 subspace. 
Intuitively, we can construct a variational wave func- 
tion by particle-hole conjugating a 6-electron Laughlin 
state with the AM — 1 edge excitation and embed- 
ding it in the IQH state with AM — 1, which wc de- 
note as m^{AM = 1)... 11010). We find that the fourth 
state in the M = 280 subspace has the largest over- 
lap (about 0.3959) with \Llj{AM = 1)... 11010); mean- 
while, its excitation energy {AE = 0.1049) is roughly the 
sum of the excitation energy of the AM = 1 edge state 
{AE = 0.07687) for the outer edge and the AM = -1 
state {AE ~ 0.02807) for the inner edge, confirming 

Eq. ca). 

Fig. [71 compares the dispersion curves and correspond- 
ing edge velocities for both the inner edge mode of the 
20-electron droplet at = 2/3 and the edge mode of 
a 6-electron Laughlin droplet at i/ = 1/3 with different 
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|(*M=45|*M=280)|^ 
|(*M=46|*M=279)|^ 
|{*i\./=47|*M=278) P 
|(*M=48|*M=277)P 
|(*M=49|*M=276)P 



1 2 3 4 I 5 

0.7527(111) 
0.6615(111) 

0.7003(1,1) 0.5814(212) 
0.7098(1,1) 0.6147(212) 0.5083(3,3) 

0.7180(1,1) 0.6238(212) 0.6000(3,3) 0.4999(5|4) 0.4130(io|5) 



TABLE II; The overlaps between the particle-hole conjugated 6-electron edge states (in 22 orbitals with d — 0.51b; Fig. [BJa)) 
and the 20-electron inner edge states (in 28 orbitals with d = 0.5/^; Fig. |6jb)). The subscript (n|m) means this is the overlap 
between the n'th state in the former subspace and the m'th state in the later subspace. 




I — ' — ^ — ' — ■ — ' — ^ ^ — ' — ' — ^ — ' — ■ — I 

-4 -3 -2 -1 1 2 3 4 



AM AM 




0.5 0.7 0.9 



d(;B) 

FIG. 7: (color online) The dispersion relation of the inner edge 
mode for 20 electrons in 28 orbitals at v = 2/3 (a) and the 
edge mode of a 6-electron Laughlin droplet in 20 orbitals at 
u = 1/3 (b) with different background confinement potentials. 
The evolutions of the edge velocities as a function of d are 
plotted in (c). It can be seen that the velocity of the counter- 
propagating edge mode in the hole Laughlin state crosses that 
of the electron Laughlin state at around d = 0.821b and the 
velocity of the electron Laughlin state has an opposite re- 
sponse to the change of the background confinement to the 
hole Laughlin state. 



background confinement. The velocity of an edge mode 
is defined asv = \de{k)/dk\. The edge excitation with an- 
gular momentum AM measured from the ground state is 
related to the edge linear momentum k = AM/R, where 
R — \/3Nel B is the radius of the A^-electron FQH droplet 
ai V — 2/3 and R = ^/QNJb for ly — 1/3. Here we smear 
the difference of the radius between the inner and outer 
edge of = 2/3. We find the velocity of the electron 
liquid Ve is larger than the corresponding velocity Vh for 
the hole droplet of the ly — 2/3 FQH state at small d 
(strong confinement), while Ve is smaller than Vh at large 
d (weak confinement). The crossing happens at around 
d — 0.821b- At d = 0.51b, the velocity of the electron 
edge is about 0.25e^/(e?i), while the velocity of the hole 



edge mode is about 0.22e^/(e?i). As d increases, Ve de- 
creases (roughly linearly) because the edge confinement 
is weaker and electrons tend to move out. However, Vh 
increases linearly with d. Therefore, the symmetry cor- 
respondence between the electron droplet and the hole 
droplet is not exact in the presence of edge confinement. 

Meanwhile, in Fig. O we pointed out that the second 
lowest eigenstate in the M = 281 momentum subspace 
is the outer edge state a.t d — 0.51b- We can expect that 
the excitation energy of the edge state decreases with 
d and may drop below that of the lowest energy state 
(presumably a bulk state) . This is indeed the case, as il- 
lustrated in Fig.m When increasing the distance d to the 
neutralizing charge background, we find an anticrossing 
behavior of the lowest two eigenstates, as their energy dif- 
ference AE = E2 — El shows a minimum at the crossing 
scale dc = 0.84/b. Overlap calculations reveal that be- 
low this dc, the lowest-energy state has a smaller overlap 
with the corresponding variational wave function of the 
edge state (jifg... 11010)) than the second lowest state. 
On the contrary, above dc, the lowest-energy state has 
a larger overlap with the variational edge state, which 
can be regarded as the ground state with an AAI = 1 
outer-edge IQH excitation. 

Fig.lSl^a) shows the edge dispersion curves for both the 
inner and the outer edge modes for different background 
confinement. Similar to the edge mode of the electron 
Laughlin state, the velocity of the outer edge mode of the 
ly = 2/3 FQH droplet decreases linearly as d increases, 
as plotted in Fig. [9jb). At d — 0.51b, the velocity of the 
outer edge mode is about 0.596e^/(e?i), slightly smaller 
than 3 times the inner edge velocity. In this case the 
small deviation from 3 suggests the two edges may be 
weakly coupled. 

We conclude the section by pointing out that we can 
identify two edge modes for a v — 2/3 FQH droplet, 
propagating along opposite directions. The outer edge 
velocity is larger than the inner edge velocity. The outer 
and inner edge modes originate from electron and hole 
droplets, respectively, and show opposite dependence on 
the strength of the edge confining potential, which breaks 
the particle-hole symmetry. 
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FIG. 8: (color online) The energy gap (red circle point line) 
AE = E2 — E\ and the overlap between the integral edge 
state |Lig... 11010) and the first two lowest states (square and 
triangular point line) in M=281 subspace as a function of 
d. The energy gap reaches its minima and the overlap has 
a crossover between the the first state and the second state 
at around dc — 0.84Z_b indicating that the first state becomes 
the outer edge state after d > dc- 



to show below, such Hamiltonians include not only the 
same hard-core interaction but also an additional one- 
body term in the electron basis. 

We start by considering a generic two-body Hamilto- 
nian in terms of hole operators: 



E 

{mi=0} 



(13) 



where the hole operators and h are related to electron 
operators: /i"*" = c and h — . It is straightforward to 
express the same Hamiltonian in terms of electron oper- 
ators: 



{mi=Q} 
- ^ Urnct^Crn + COnSt., 



where 



{Vr 



nrnnm 



Vr 



nmmn 



} 



(14) 
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FIG. 9: (color online) The dispersion curves of a = 2/3 FQH 
droplet with 20 electrons in 28 orbitals (a) for both the inner 
and outer edge modes at different d's. As in the 6 electron 
Laughlin state, the velocities of the outer edge mode decrease 
linearly as a function of d and have larger velocities than the 
inner edge (b) in the whole parameter range where the global 
ground state resides in the M=280 subspace. 



VI. PARTICLE-HOLE TRANSFORMATION 

The Laughlin state is the exact zero-energy state of 
a special two-body Hamiltonian with hard-core interac- 
tion. In this section we use particle-hole transformation 
to construct Hamiltonians that make the hole Laugh- 
lin states (which we used as variational ground states in 
previous sections) exact ground states. As we are going 



is the Hartree-Fock self-energy of the state m when all 
the Norb orbitals are occupied by electrons. Thus, in the 
electron basis we get the same two-body interaction plus 
a one body potential, which is attractive if the two-body 
potential is repulsive. 

From now on we focus on the specific short-range 
(hard-core) interaction that corresponds to Haldane 
pseudopotential Vm — <^i,m for Hh- After diagonalizing 
the Hamiltonian Hh with TV = 20 electrons in Norb — 28, 
we obtain the energy spectrum in Fig. 1101 It is worth 
pointing out that since there is no edge confinement other 
than the momentum cut-off due to the choice of Norb, we 
have Nh = Norb — N = 8 holes in the system. Not surpris- 
ingly, the largest angular momentum of the degenerate 
ground states is Mq = 294, as expected from Eq. 
with Ni = Norb = 28 and A^e = 20. This is the densest 
ground state configuration for holes and is thus incom- 
pressible. For AM = M~Mo = -1, -2, -5, we find the 
ground state degeneracy in each subspace is n{AM) = 1, 
2, 3, 5, and 7, respectively. This is precisely the number 
expected for the Laughlin droplet (except for this case, 
the momentum is negative), as generated by 



00 ^ 

5:.(AM),-.nT3V 

A 7l ,f ™ — 1 y 



(16) 



AM 



In the spirit of Refs. 33 and 34, a mixed Hamiltonian 
that contains both the Coulomb Hamiltonian He used 
in earlier sections and the above hard-core Hamiltonian 
Hfi parameterized by A is considered: 



H = XHh + (1 - \)H 



c- 



(17) 



The idea here is that the hard-core Hamiltonian H^ raises 
the energy of bulk excitations while having little effect on 
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FIG. 10: (color online) The energy spectrum of the P-H con- 
jugated Hamiltonian with hard-core interaction (Eq. (|14|) ') for 
20 electrons in the 28 orbitals. M is the total angular mo- 
mentum for the 20 electrons. The first 50 energy levels are 
plotted for each M. The lowest energy states (they should be 
the exact zero energy state if we include the constant term) 
at M=294,...,289 are degenerate with degeneracy 1,1,2,3,5,7. 

the edges states, thus its presence helps separate the two 
energetically. Fig. [TT] shows the energy spectrum for 20 
electrons in 28 orbitals in the pure Coulomb case (A = 0) 
and a mixed case (A — 0.5). Although the ground state 
of Hh for 20 electrons in 28 orbitals is the 8-hole Laugh- 
lin state, for the mixed Hamiltonian above we obtain a 
different ground state with the same quantum number as 
the 6-hole Laughlin state {M = 280). This is because the 
neutralizing background charge at d = 0.9 serves as a re- 
pulsive potential to the holes that pushes two holes to the 
outer edge. It is clear that edge excitations show up at 
lower energies for the case A = 0.5. In particular, overlap 
calculations indicate that in the case of pure Coulomb in- 
teraction with A = 0, the third state in M=280 subspace 
is the simplest combination state which has the largest 
overlap (0.194) with \L\,j{/^M = 1)... 11010), while in 
the case of the mixed Hamiltonian with A = 0.5, it is 
the second state in M=280 subspace that has the largest 
overlap (0.296) with \L\,j{AM = 1)... 11010). Therefore, 
the mixing of the hard-core Hamiltonian indeed has the 
effect of separating the band of edge states from the bulk 
states. 



FIG. 11: (color online) The energy spectrum for the 
pure Coulomb Hamiltonian with A = (a) and a 
mixed Hamiltonian with A = 0.5 (b). When A = 
the overlaps between the outer edge states (and ground 
state) and their corresponding conjugated states are: 

|(*M=280i|Lf6>il' = 0.684, K*M=281i|if6-11010)P = 

0.385, K*Af=2823|Lf6--10110>|^ = 0.133, 

i(*M=282i5|i?6-11001)p = 0.158; the over- 
lap for the simplest linear combination state is: 
i(*Af=2803|i?7(AM = l)...11010)p = 0.194. In the 
case of a mixed Hamiltonian with A = 0.5, they are: 

K*M=280i|Lf6>igP = 0.66, |(*Af=281i|if6-11010>P = 

0.539, K*M=2822|i!6-10110>p = 0.218, 

l(*Af=282io|i?6-11001)|" = 0.107, |(*M=2802|i?7(AM = 

1).. .11010)1^ = 0.296. 




FIG. 12: (color online) The dispersion curves for both the 
inner edge and the outer edge mode for 20 electrons in 28 or- 
bitals with d — 0.91b and different layer thicknesses. The the 
layer thickness softens the interaction between the electrons 
and reduces the edge mode velocities. 



VII. LAYER THICKNESS 

One improvement in the semi-realistic model is to con- 
sider the effect of finite electron layer thickness. In real 
experimental samples, quasi-two-dimensional electrons 
are confined in the GaAs quantum well, whose width can 
be as large as 30 nm or a few magnetic lengths. Since 
the vertical motion of electrons is suppressed at low tem- 
peratures, the quasi-2DEG can be approximated, to the 
lowest order, by an ideal 2DEG located at the peak of the 
wave function in the perpendicular direction. The finite 
width softens the repulsion between electrons and thus, 
together with other factors (like higher Landau level). 



may help stabilize certain fragile FQH states . ' In this 
section, we bricfiy discuss the effects of the 2DEG layer 
thickness on the velocities of = 2/3 edge modes. We 
use the Fang- Howard'*'^'^"'^ variational wave function, 

Zoiz) = 2(26)-3/2ze-"/2fc^ ^^g) 

to model the electron layer thickness effect, where 6 is a 
measure of the well thickness. We obtain the same qual- 
itative behavior for an infinite quantum well potential. 

We can integrate the Fang-Howard wave function to 
obtain the renormalized Coulomb interaction in Fourier 
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space 



VFH{k) 



1 3(fcfe)2 + 9fc6 + 8 
8fc {kb+lf ■ 



(19) 



In Fig. [121 we show the edge dispersion curves for both 
the inner and the outer edge modes for a = 2/3 FQH 
droplet of 20 electrons with different layer thicknesses for 
d = 0.9. The ground state angular momentum is again 
M = 280, consistent with the 6-hole Laughlin droplet 
picture. We find that the velocities of both edge modes 
are reduced by increasing layer thickness, as expected. 
This is in contrast to the effects of the confining potential, 
which modifies edge velocities in an opposite way. 



VIII. QUASIPARTICLE 

In this section, we demonstrate that both one quasi- 
hole and one quasiparticle can be excited at the center of 
the FQH droplet with an additional short-range impurity 
potential. This section is a natural generalization of sim- 
ilar works by some of the authors for the Laughlin case 
aX u — 1/3^ and for the Moore-Read Pfaffian casci^i^ 

We follow the previous work^S by using a Gaussian im- 
purity potential Hw = Wg J^m 6xp(— m^/2s^)c^Cm with 
a finite width s = 21b to excite and trap either a quasi- 
hole or a quasiparticle. We consider a system of 20 elec- 
trons at = 2/3, whose ground state angular momentum 
is M = 280 for d = O.bh- At Wg = 0.2 and Wg = -0.2, 
the global ground state resides in the angular momentum 
subspaces at M = 286 and M = 274, respectively. The 
change of the ground state angular momentum M by ±6 
suggests that we have induced a charge |e|/3 quasihole or 
a charge — |e|/3 quasiparticle at the center of the 6-hole 
droplet. The density plot for the electron ground state in 
Fig. [T3] confirms the charge depletion and accumulation 
in the corresponding cases. The increase of the electron 
ground state angular momentum means a corresponding 
decrease in the angular momentum of the hole droplet, 
suggesting a quasiparticle excitation for the hole droplet 
and thus a quasihole excitation for the v = 2/3 electron 
ground state, as plotted in Figs. \W( c) and (d). The de- 
crease in angular momentum, on the other hand, suggests 
a quasiparticle excitation, as plotted in Figs. fTSTe) and 
(f). These excitations are localized at the origin, where 
we apply the impurity potential, and their presence has 
no effect on the edge excitation spectrum, since these are 
Abelian anyons. Roughly speaking, the quasihole and 
quasiparticle excitations induce the same density pertur- 
bation on the ground state, except for the opposite signs, 
suggesting a quasiparticle-quasiholc symmetry. 



IX. CONCLUDING REMARKS 

To summarize, we study the ground states, edge, and 
bulk excitations of = 2/3 fractional quantum Hall con- 
densates in a semi-realistic microscopic model. We find 




FIG. 13: (color online) The density profile of the 2/3 state 
(1/3 hole-Laughlin state) and its quasihole and quasiparti- 
cle excitations. The system contains 20 electrons in 30 or- 
bitals with d — 0.51b and the Gaussian tip potential Hw = 
Wg exp(— m^/2s^)CmCm has a width s — 21b- (a) and (b) 
are the density profiles for the ground state with angular mo- 
mentum M=280, (c) and (d) are the density profiles for one 
quasihole state with M=286 when Wg = 0.2 and (e) and (f) 
are for one quasiparticle state with M=274 and Wg = —0.2. 



strong numerical evidence that a v = 2/3 droplet can 
be regarded as a i^/i = 1/3 Laughlin hole droplet embed- 
ded in a larger lyj — I integer quantum Hall droplet. In 
particular, we find two counter-propagating edge modes, 
which are associated with the inner edge (the edge of the 
hole droplet) and the outer edge (the edge of the integer 
quantum Hall edge). The inner edge mode is well sepa- 
rated from bulk excitations and resembles the edge of a 
1/ = 1/3 electron droplet, except that they propogate in 
opposite directions and respond oppositely to the edge 
confining potential, which explicitly breaks the particle- 
hole symmetry. The outer edge states have higher energy 
and mix with bulk excitations due to the computational 
limit on the Hilbert space dimension. The v = 2/3 quan- 
tum Hall droplet also has the same ±e/3 quasiparticle 
and quasihole excitation as a, v — 1/3 Laughlin droplet. 
These features are robust in the presence of finite electron 
layer thickness, which softens the Coulomb interaction 
between electrons. 

One of the major advantages of the disk geometry is 
that we can identify edge modes and determine the veloc- 
ities of the edge modes. We have previously applied the 
same method to — 5/2 fractional quantum Hall systems 
and found significant differences in charge and neutral 
velocities;^ which leads to quite different temperature 
regimes in which charge e/4 and charge e/2 quasiparticles 
can be observed in interference experimentSf2Lil Simi- 
larly, based on comparison with the edge excitations in 
the ^ = 1/3 Laughlin case and in the v = I integer case. 
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we are able to resolve the edge excitations of the = 2/3 
case. It is interesting to point out that the outer edge 
velocity is roughly 3 times that of the inner edge mode 
velocity, which is about the same as the edge mode veloc- 
ity of the V — l/i Laughlin liquid with the same Coulomb 
interaction and a similar confining potential strength due 
to the neutralizing charge background (Fig. [3:). In gen- 
eral, edge mode velocities are non-universal, depending 
on details of electron-electron interaction and confining 
potential. For v = 2/3, at which there are two counter- 
propagating edge modes, the velocities further depend 
on the coupling between the two edge modes affected by 
interactions and impurities j^"^'^^ for example. The fact 
that we are observing a relatively robust ratio (~3) of 
the outer edge mode velocity to the inner one (consistent 
with their density changes) suggests that, at the length 
scale of our finite-size study, the velocities are dominated 
by the Coulomb interaction strength determined by the 
electron density change associated with each mode and 
that the two edge modes are very weakly coupled. We 
note that, in the thermodynamic limit, the long-range 
nature of the Coulomb interaction is expected to force 
the two modes to reorganize into a charge mode and a 
neutral mode, with the charge mode velocity logarithmi- 
cally divergent while the neutral mode velocity remains 
finite in the long wave length limit. The systems size of 
our study is too small to see this trend. 

Another feature of the current calculation is that we 
consider a semi-realistic confining potential arising from 
the neutralizing background charge. This, together with 
long-range interaction between electrons, allows us to 
compare the energetics of competing states and discuss 
the qualitative dependence of eigenenergies and edge 
mode velocities on the confining potential. This is an 
extremely interesting subject, especially when we have 
the = 5/2 quantum Hall systems in mind. In that 



case, there are at least two competing candidates for 
the ground states, the Moore-Read Pfaffian state and its 
particle-hole conjugate, the anti-Pfaffian state. They are 
exactly degenerate, if we neglect the particle-hole sym- 
metry breaking terms, such as the 3-body interaction due 
to Landau level mixing. As demonstrated in this study, 
the confining potential is also a relevant symmetry break- 
ing term, which leads to an opposite dependence of the 
edge mode velocities on this potential. It has been pre- 
dicted'^'' that the anti-Pfaffian state is favored in weak 
confinement (smooth edge) while the Moore-Read Pfaf- 
fian state is favored in strong confinement (sharp edge). 
It is also interesting to point out that for v — 2/2> the de- 
pendence of the edge mode velocity of the inner edge on 
the confining potential is significantly weaker than that 
of the outer edge (Fig. [5]), suggesting a screening effect 
by the outer edge on the inner edge to the change of the 
confining potential. 
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